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Abstract 

> 

^ ■ We consider the de Sitter cosmology deformed by the presence of a thermal bath 

^ . of radiation and/or time-dependent moduli fields. Depending on the parameters, 

^ i either a first or second order phase transition can occur. 

^ ■ In the first case, an instanton allows a double analytic continuation. It induces 

. a probability to enter the inflationary evolution by tunnel effect from another 

^ ' cosmological solution. The latter starts with a big bang and, in the case the 

transition does not occur, ends with a big crunch. A temperature duality ex- 
^ ■ changes the two cosmological branches. In the limit where the pure de Sitter 

' universe is recovered, the tunnel effect reduces to a "creation from nothing" , due 

to the vanishing of the big bang branch. However, the latter may be viable in 
some range of the deformation parameter. In the second case, there is a smooth 
evolution from a big bang to the inflationary phase. 
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1 Introduction 



The recent astrophysical observations indicate that our universe is in a phase of classical 
expansion, with a small positive cosmological term, A, representing 60 — 70% of the total en- 
ergy density. However, at early time, quantum corrections to this trajectory are expected to 
be of first significance. For instance, drastic non perturbative effects could occur as topology 
changes, as can be seen in the context of field theory in semi classical approximation. As 
an example, a topology change is provided by the instability of the five dimensional Kaluza- 
Klein (KK) space-time [1]. In that case, one observes the nucleation of a finite size "bubble 
of nothing" that is then growing up at the speed of light. The transition is described by tun- 
nel effect in terms of an instanton configuration allowing a "double analytic continuation". 
This means that at two different Euclidean times Xj and t/, analytic continuations to real 
times T = Ti + it and r = tj + it are allowed. To be specific, the configuration admits a time 
independent asymptotic behavior allowing a continuation at = +oo to the KK universe, 
while another continuation at tj = describes the evolution of a bubble. 

Actually, the KK universe is suffering from a first order phase transition where bubbles 
appear instantaneously, grow and coalesce, so that the space "evaporates" into nothing (after 
an infinite time, since the volume is itself infinite). In some sense, some reversed ideas can 
be invoked in another example of topology change, namely in Vilenkin's scenario of the de 
Sitter space "creation from nothing" [2,3]. In that case, a finite radius space appears 
instanteously by tunnel effect from a space-time state that amounts to the empty set. This 
5*^ "bubble" is then following a de Sitter growing up evolution. The transition involves 
an instanton, whose shape is an S*^ hemisphere. An analytic continuation on its boundary 
amounts to "gluing" a de Sitter universe, while the fact that the instanton is compact with 
no other boundary to analytically continue allows an interpretation in terms of a transition 
from an empty set. 

Concerning the birth of the de Sitter space, the instanton method provides an estimate 
of the transition probability equivalent to the one derived from the Hartle-Hawking wave 
function [4] approach. The transition amplitude and wave function being proportional to 
e""^^, where Se is the Euclidean action, the probability of the event is 

p oc e~2^^ . (1.1) 
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A selection principle of the cosmological constant A based on the extremum of p leads to a 
favored value A — > 0+, [5]. However, a present too small cosmological term cannot account 
for the 60% to 70% of the total energy density, which is necessary to explain the recent 
dark energy data. To remedy to this problem, it as been stressed [6, 7] that since the de 
Sitter space has an horizon, it can be associated a Hawking temperature. Thus, quantum 
fluctuations of the metric (or any massless mode introduced in the model) induce a space 
filling thermal bath. The latter implies an additional radiation term in the action^ and thus a 
back reaction on the metric background that deforms the de Sitter solution (see also [9]). In 
that case, the modified probability transition derived from the Euclidean action is maximal 
for a non vanishing finite value of A. This computation has also been addressed by the 
authors of [10] and refined in [11], who considered the Wheeler-de Witt equation in the WKB 
approximation. They found that the tunnel effect is not connecting the deformed de Sitter 
expansion to "nothing" , but to what they called a thermally excited era. However, the latter 
is A-dependent. This point makes the difference with the pure de Sitter space created from 
a A-independent state, namely the empty set, so that extremizing the transition amplitude 
with respect to A was making sense. In the case of the deformed solution, applying this 
selection rule for A is thus questionable. 

Our first aim is to reconsider the thermally deformed de Sitter space-time. The previ- 
ous analysis that involves an instanton transition occurs when the radiation energy density 
is below some critical value. This case corresponds to a first order phase transition. A 
double analytic continuation on the Euclidean solution implies a probability to connect two 
cosmological behaviors in real time. The first one is the de Sitter like evolution, while the 
second one is an era that begins with a big bang and ends with a big crunch, when the 
transition does not occur. If instead the radiation energy density is above the critical value, 
the phase transition becomes second order. It describes a smooth evolution from a big bang 
to a de Sitter like behavior. In all these solutions, the temperature is formally infinite when 
the radius of the universe vanishes. This means that these evolutions should be trusted as 
soon as the temperature passes below some upper bound such as the string or Hagedorn 
temperature Th of an underlying fundamental string theory. 

We then generalize the model by including time dependent moduli-like fields that im- 

"'^See [8] for a cosmological scenario based on a cascade of transitions between vacuum and radiation 
energy. 
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ply a further deformation of the de Sitter solution. Since, these deformations respect the 
Robertson- Walker (RW) isometries, they can be expressed in terms of the scale factor a of 
the universe. The Hubble parameter is taking the form 



where the constant and l/o? terms are associated to the cosmological constant and uniform 
space curvature. The monomial accounts for the effective radiation term discussed 
before, while the on shell moduli contribution reaches the 1/a^ contribution. What we find 
is that the moduli deformation can be fully absorbed in a redefinition of time. We then 
describe different cosmological scenarios that depend on Cr and Cm- These parameters can 
be expressed in terms of the cosmological constant, the number of effective "massless" degrees 
of freedom, the temperature at the transition and the slope of the moduli dependence. If this 
two parameters model is interesting by itself from a formal point of view, its study is also 
motivated by the fact that it describes more complicated systems once they are considered 
on shell. For instance, cases involving a scalar field with a non trivial potential and coupled 
to moduli with non trivial kinetic terms are treated in [12]. 

2 Effective cosmological models in minisuperspace 

We focus on models where space is 3-dimensionaP and closed. The number of degrees of 
freedom is important, and actually infinite in the context of string or M-theory. However, a 
subset of them to be considered in an effective field theory is determined by the ultraviolet 
cut off set by the fundamental theory. Some of the simplest class of models that can be 
constructed are considering the so called "minisuperspace" (MSS) approximation [4, 13]. 
They involve isotropic and homogeneous spaces only, so that the only gravitationnal degree 
of freedom is the scale factor of the universe, to which some matter fields can be coupled to. 
The space-time metric is thus 



^By this, we mean that if the fundamental theory is for instance string (or M-)theory, there are compact 
dimensions whose sizes are dynamically constrained to remain of order Is , the string length. 





ds'^ = {-N^{t)dt^ + a^{t)dnfj , where 
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Mp is the Planck mass and aa{t) is the radius of a 3-sphere. is a gauge dependent function 
that can be arbitrarily chosen by a redefinition of time'^. (In our conventions, and t are 
dimensionless.) In order to include in the effective action the quantum fluctuations of the 
full metric and matter degrees of freedom of the theory, we switch on a thermal bath which 
is created simultaneously with the space-time. The thermal bath involves the degrees of 
freedom below the temperature scale T and is consisting in an extra radiation term in the 
effective MSS action. This amounts to take into account the back reaction of the thermal 
bath on the space-metric i.e. deform the de Sitter solution [6,7]. 

The starting point is a gravity action coupled to some matter fields, whose Lagrangian 
density is Lm '■ 

S = Sg + S^ = ^J d\^g{R - 2A) + y d'x,/^L^ . (2.2) 

Assuming the metric (2.1) in the gravity part and freezing all other degrees of freedom, we 
obtain the MSS-action, 




a a 



N 



, where A = —a . (2.3) 



In this expression, a and A^ have been chosen positive, without loss of generality. Further- 
more, the boundary terms arise from the integration by parts of the second derivative of a 
in the scalar curvature R. 

The non-vanishing components of the energy-momentum tensor which are induced by 
the matter fields (and which are invariant under the isometries of the RW space-time), can 
be expressed in terms of the energy density p and pressure P : 

Tu = pa^N' , Tij = P<J^a''~g,, , (2.4) 

where drt\ = Qijdx^dx^ is the line element on S^. These quantities appear in the equations 
of motion of A^ and a i.e. the time-time and space-space Einstein equations : 

+ 3^ = Uir'a^N^p + 3XN\ (2.5) 



^We shall use in the following the time variable t either when the function N{t) remains unspecified, 
or when it is fixed to N{t) = 1. We shall instead use another variable, i, when we choose to fix N{i) ~ 



1/ (l + A ^K/a^^i)^, where k is a constant. 
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a a N f N'^ 949 ,9 , ^ 

-2- + 2-— - ^ = Un^a'^N^P - 3\N^ . 2.6 

a a i\ \a J a'^ 

Since can be fixed to any arbitrary function by a reparameterization of t in tlie metric 
(2.1), this differential system of two unknown functions must be degenarate so that we only 
have to solve any of the two equations^, e.g. the Friedman one (2.5). Let us consider some 
cosmological solutions associated to various matter contents and where phase transitions can 
occur. 



3 Pure de Sitter space 

With neither matter nor thermal bath, we have p = P = 0. Choosing N{t) = 1, eq. (2.5) is 
easily solved, 

a{t) = (3.1) 

This corresponds to an usual de Sitter evolution where a phase of contraction from an 
infinitely past time is followed by an expansion phase when the radius a has reached its 
minimum. 

However, a different cosmological scenario can give rise to an identical de Sitter expansion. 
Actually, one can consider the real time evolution for t > as the result of a tunnel effect, 
since an analytic continuation at tj = —irf = is allowed. The instanton solution in question 
is obtained by substituting t = —ir in (3.1)^, 

|cos(yAr)| 

= ' (3.2) 

where we find convenient to insist on a conventional positivity of a^; for future purpose. The 
Euclidean metric takes the form ds'^ = ((T^/A)(Ac/r^ + cos'^{\/XT)dQ.l), so that the instanton 
is a 4-sphere of radius a/V^, if \/Xr G [— 7r/2, 7r/2]. 

Vilenkin, assuming that the range of the Euclidean time [rj, rj] = [— 7r/2-\/A, 0], considers 
this instantonic hemisphere as giving rise to a transition amplitude from nothing (a state 
of the space-time where it is the empty set) to an space with de Sitter cosmological 
evolution [2,3]. As we shall see later, this interpretation is very natural in the deformed de 



"'This is not true for static solutions. 

^For any field /(t), we define in general its Euclidean counterpart fsir) = /(—ir) and its derivative 
fsiT) = -ifi-ir). 
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Sitter solution which is induced by the thermal bath. The Euclidean action, 



a\aE 



(3.3) 



evaluated on shell, has vanishing boundary terms and takes the value 

SEg = -^. (3.4) 

so that the probability transition 

p oc e2/(3A) (3.5) 

becomes maximal for a vanishing cosmological constant. Considering the second hemisphere 
of the S'^-instanton, \/~\Ti = and VXtj = 7r/2, we note that the previous probability also 
evaluates the chances for a contracting de Sitter space to simply vanish instantaneously when 
it reaches its minimal radius at t = 0. 

Another derivation 

Before closing this section, we would like to mention an alternative method to solve the 
Friedman equation (2.5) that will appear useful to express the general deformed de Sitter 
solution once the thermal and moduli effects are simultaneously present. Consider 

N'ii) = ^ ^ J . , , (3.6) 
1 + A ^K,la'^[t) 

and choose k = — 1 so that the 1/a? term of eq. (2.5) is cancelled out^, 

*)^A. (3.7) 

a J 

There are two solutions to this equation, whose corresponding metrics are 

defined for t > — log(A)/2V^ and i < log(A)/2\/A, respectively. These parameterizations 
do not admit analytic continuations to Euclidean time, and actually correspond to "half 
solutions" , as can be seen as follows. Integrating in each case 

dt = N{i)di, (3.9) 



^As signaled in footnote 3, remind that the "dot" derivative is here with respect to t, the time variable 
in the gauge (3.6). 
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one can choose the time variable t (corresponding to the gauge = 1) such that 



(3.10) 



to glue the two branches in t in a smooth way, t G M (see fig. 1). In that case, the metrics 
ds"^ are taking a RW form in the = 1 gauge (see eq. (2.1)), with a defined in eq. (3.1). 



Fi gure 1: As a function oft, the pure de Sitter scale factor a has two branches. As a function oft, these 
branches are glued together in a smooth way. 

4 Thermally deformed de Sitter universe 

The quantum fiuctuations of the previous background metric have been shown to induce 
an additional effective "radiation term" in the action. This implies a back reaction of the 
massless gravitationnal degrees of freedom on on the metric [6,7]. In other words, the solution 
of the new equations of motion describes the evolution of a RW closed universe in presence 
of thermalized radiation. An explicit derivation of the effective thermal contribution can 
be found in [6,7], while the consistency of the semi-classical approach has given a similar 
result [10]. We would like here to reconsider these effects by using thermodynamical and 
precise dynamical constraints to reveal new features. The results can be applied to any 




a(f) 



am 
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model involving massless bosonic or fermionic modes (beside the gravitationnal ones), as 
long as their vaccum expectation values do not show up in the effective action. 

For bosonic (or fermionic) fluctuating states of masses nib (or nif) in thermal equilibrium 
at temperature T, the general expressions of the energy density pt and pressure Pt are 

pT= ip(^^)+ E ^pi^^fh pt= E ^^(^/)' (4-1) 

boson b fcrmion / boson b fermion / 

where 

and E{mb(^f)) = + ^1^^-^ . If the system is consisting in nB[F) massless degrees of freedom, 
it is easy to see that the state equation 

Pt = ^Pt = ^ {nB + (4.3) 

is satisfied. These relations combined with the expression of the energy-momentum conser- 
vation 

PT + 3 {pT + Pr) = (4.4) 

is then implying that 

aT = constant (4-5) 
^2 I 

PT = , „ 2 4 \ — ' (4-6) 

where 5|,/(167r^cr^A) is a positive integration constant. The factor of A in its denominator is 
chosen for later convenience. 

To proceed, we solve the Einstein equations in the case where the energy density and 
pressure are the thermal ones. This has been done for eq. (2.6) in [6,7,10] and [14] in the 
N{t) = 1 gauge. However, we prefer to deal with eq. (2.5) as in [9] but in a way similar to 
what we did at the end of the last section for the pure de Sitter case, since this will be easier 
to generalize when we take into account moduli fields. For this purpose, we choose again 
N{t) with a form 

= } , (4.7) 

1 + A ^Kia^[t) 



and that 
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and look for k so that the thermal contribution is cancelled out. This implies that 

+ K + ^ = , (4.8) 

a condition that admits real solutions as long as the constraint < 1 is satisfied. 
The case 5'^ <1 

The Friedman equation (2.5) becomes 

a\ ^ K + 1 , , , 

- +^ = A, 4.9 

showing that the solution of our present problem can be related to the one of section 3, the de 
Sitter case with no thermal effect, up to a shift in the curvature term. The hyperbolic 
cosine solutions of eq. (4.9) for the two roots of (4.8) are giving rise to the metrics 



dsl = ( + ^^'{^ cosh\v^i)dnl] , (4.10) 

\ cosh2(v^t) / 

where 

u= ]-'^^ e]0,i], (4.11) 

if > 0. Since they simply correspond to different gauge choices of functions A^, we 
already know that they are equal up to a redifinition of time. This means that the potential 
singularity of ds"^ when the denominator of N(t) vanishes should be fake. In any case, let 
us make contact with the more intuitive = 1 gauge by choosing a new time variable t 
satisfying eq. (3.9). This will in particular allow us to study in an explicit way various 
analytic continuations. 

For ds"^, one can fix the integration constant such that 

cosh2(v^t) = n + (1 - n) cosh^i^/Xt) , (4.12) 
so that the metric is taking the form ds^ = (— (it^ + {i(t)^ dfl'^), with 



a = Ar^e + cosh^{VXt), (4.13) 

where 



For ds'^, one can have instead 

cosh2(v/At) = 1 + - 1) cosh^{\^t) , (4.15) 

which has the effect to glue the two branches -y/At > arccosh(u~^/^) and ^/Xi < — arccosh(n~^/^) 
in a smooth way, with t G M (see fig. 2). This time variable, has the property to explicitly 
show that ds"^ = ds\. 




Figure 2: When one is considering the time variable t involved in ds^_ in the pure thermal case, the scale 
factor a has two branches. In terms of t, these branches arc glued together in a smooth way. At the minimum, 
one has a — a+ given in cq. (4.18). 



This solution is a deformation of the de Sitter solution (3.1) i.e. a contraction phase 
followed by an expansion one. As before, the latter can also arise by tunnel effect since 
the analytic continuation at ti = —iTf = is still valid. Let us consider for a moment the 
relevant instanton solution obtained under the substitution t = —ir, 

ttE = J\f\J e + cos^ (v^r) , (4.16) 

for the range of Euclidean time — 7r/2-\/A < r < 0. An important point is that at ^/Xt^ = 
— 7r/2, a^; is no longer vanishing as in the birth of the de Sitter space (ref. [15] presents 
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another example of this phenomenon). Instead, a second analytic continuation is allowed, 
\/\t = —7t/2 + i\/\t, and another cosmological solution 



a = U\Je-smh^{\/\t), (4.17) 

satisfies (2.5) (in the = 1 gauge), for — arcsinhy^ < \f\t < 0. Actually, this sinh-solution 
is connected to the cosh-solution (4.13) via the instantonic one (4.16). The cosmological 
evolution makes use of all above mention solutions. 

• The universe starts from a big-bang era at t ~ — A^^^^arcsinhA/e, with an evolution fol- 
lowing the sinh-solution till t = 0. 

• At t = 0, either i) the space starts to contract till a big crunch occurs att = A~^/^arcsinhy^, 
or a) there is a first order phase transition that changes instantaneously its radius. 

• In case ii), the cosmological evolution for later times, t >0, follows the infiationary cosh- 
solution. 

Actually, the very early part of these scenarios (and late time in case i)) are not trustable, 
assuming that the fundamental theory is string like. The degrees of freedom grow exponen- 
tially at high temperature, and blow up at the Hagedorn temperature. This effect gives a 
cut in time before which the sinh-solution is no more valid (and another cut for late time in 
case i)). A stringy approach is then necessary to cover such ultra hot eras of the universe. 

At the transition, the scale factor is jumping. 



1 -0^ 
2A 




^ a+ = \ r ' \: , (4.18) 



(see fig. 3). The transition probability at t = (see section 5.2) depends on A and on the 
deformation parameter e (or equivalently on 6^). The maximal value of this probability is 
for finite A and depends on e. This is fundamentally different from the pure de Sitter case 
where the maximal probability is for A = 0. 

It is interesting to make a comparison between the pure de Sitter case and the thermally 
deformed one. When 6^ 0, the big bang occurs at tj ~ —6t/'2^/X, while we also have a_ ~ 
5t/2v^ at t = 0_. This means that Vilenkin's interpretation of the spontaneous creation 
of the pure de Sitter space can be seen as the limit when the big bang sinh-cosmological 
branch has shrunk to nothing (see fig. 4). Actually, the instantonic jump a_ a+ becomes 
more drastic as the thermal effects become more negligible, 5t ^ 1, since a+ ~ 1/a/A in 
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Fi gure 3; Cosmological scenario of an expanding RW universe in presence of thermal and time- dependent 
moduli effects below some critical value. A phase of expansion starts from a big bang and is connected to an 
inflationary deformed de Sitter one by an instanton. This describes a first order transition, where the scale 
factor jumps from a_ to a+ . 



this limit. On the contrary, the transition from the big bang branch to the cosh-inflationary 
deformed de Sitter one is smoother and smoother as we approach the critical magnitude of 
the thermal effects, 5t ^ 1. Actually the amplitude of the instanton solution decreases as 
we approach this value of the parameter. When 5^ = 1, the instanton solution is constant, 
implying that its analytic continuations are static, a{t) = l/v^2A. 

An estimate of the transition probability from the big bang branch to the deformed de 
Sitter one is p oc e"^"^^"-**", where 

s.., = 1 £' iv.4 (~ (^) ' - 4 + A + Jj) + 1 

is the Euclidean effective action that is giving rise to the desired equations of motion. In 
this expression, SEeff is evaluated with the solution (4.16) and has vanishing boundary terms 
for -tt/2 < V\t < 0, (see eq. (5.57)). However, one could sum in principle the instanton 
contributions for the series of ranges —{2n+ l)7r/2 < ^/Xt < 0, n G N. The total probability 
would take the form p^^j = X]n>o '^sn+ie"'-^"^"'^''^'^-^"-*'', where C2n+i are positive constants. Sim- 



(4.19) 



12 




-TZ/lX 



111 



Fi gure 4: Thermal instanton (in hold) compared to the pure de Sitter S one. In the limit where thermal 
effetcs are negligible, 5t "C 1, the instantonic jump of a is drastic, passing from a_ ^ St/2^/\ to 0+ ^ l/\/X. 
When thermal effects are getting closer and closer to their critical magnitude, 6t ~ the amplitude of the 
transition from the big bang branch to the inflationary deformed de Sitter one is decreasing. At the very 
precise value 5t — 1, the instanton solution is a constant, aE{T) = 1/V2X. 



ilarly, the probabihty to pass from the contracting phase of the deformed de Sitter space to 
the expanding one would take the form p[^^ = X]n>o ''2ne~^"'^^*=-^, for positive constants C2n, 
(see fig. 5 and [16]). 

Introducing the temperature at the transition 

It is useful to express the deformation parameter 5|> (or Af and e in eq. (4.14)) in terms 
of more intuitive quantities, such as temperatures. First, observe that before and after the 
transition, the quantity aT remains constant (see eq. (4.5)). Since these evolutions are 
related by a double analytic continuation, the constants must be equal. This allows us to 
express aT in terms of either the temperatures T+ at t = 0+ or r_ at t = 0_ : 



1 ± 



1-52 , 



(4.20) 



Plugging either of these relations in the expression of px in eq. (4.3) and comparing the 
result with eq. (4.6) imposes a consistency condition to be satisfied by 6^. Its solution 
involves a mean temperature Tm satisfying T+ < Tm < T_, 

1/4 



6^ 



iTi/n+n/Tir 



where T„ 



15 



A 



Inp) 



1 
a 



(4.21) 
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Fi gure 5: The inflationary phase of the deformed de Sitter space for t > can arise from any contracting 
phase or big bang phase for t < 0, after a multi-instantonic transition. 



This implies that T^T- = and that 5^ is invariant under the duahty transformation 

(4.22) 

Actually, this operation exchanges the scale factor of the thermally deformed de Sitter evo- 
lution, 

^ l + ^^cosh(2yAt) , t>0, (4.23) 



with the big bang / big crunch one 



a 



,,1+ 7 - cosh(2v^t), t<0. (4.24) 
2A y "T - 

In some sense, the instantonic transition between them is "inversing" the temperature, since 
T(t) > T_ for t < 0, while T{t) < T+ = T^/T_ for t > 0. If one wishes, one can also rewrite 
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the constants M and e of eq. (4. 14) in terms of either T+ or T_ 




Af=—, r e = ^ (4 25) 

and note that the function A^^ appearing in eq. (4.7) satisfies N^{t) = 1/(1 — T'^{t)/T^) for 
ds"^, respectively. 

Focussing on the X-dependance of 6^ 

We would like here to consider the dependance on A of This is relevant, for instance, if 
one wants to minimize the Euclidean action with respect to A (keeping fixed the number of 
degrees of freedom and temperatures T±). Doing so was considered as a selection principle 
for the cosmological constant in [6, 7, 10, 11, 17]. The reason for this was based on the 
fact that this action (see eq. (5.57)) is controlling the probability transition between the 
two cosmological evolutions. However, there are two hypothesis for this argument to be 
considered. The first one is that the viable universes are the deformed de Sitter ones, i.e. 
that the big bang / big crunch evolution is too hot or too short in time. However, this 
condition implies that 6^ is small enough, which in turn can be translated into a condition 
on A itself (see eq. (4.21)). The second is that the probability to start the A-dependent big 
bang evolution is almost constant clS db function of A. 

Eq. (4.21) can be rewritten as 

4 = ^ ^—.=- (4.26) 

A (1 + jy+/\f (1 + Xjv^f 

where 

(nB + -np)Tla\ (4.27) 
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and reproduces the result of [11]. Actually, A is the geometric mean of z/+ and 

so that eq. (4.26) is relevant to derive an expansion of 6^ for small t^+fX, whose leading term 
reproduces the result of [6,7,10]. Eq. (4.26) goes beyond the small z/+/A approximation by 
taking consistently into account the back reaction term 5|, that appears in the r.h.s. of eq. 
(4.20). 
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A conceptual difference of our approach compared to that of [6, 7, 10] is that we insist 
on the difference between the two temperatures T±, where T_ > T+. Thus, the only stringy 
constraint to not reach the Hagedorn temperature during the cosmological evolution is T„ ^ 
Th- In that case, the temperature T{t) remains small, T{t) ^ T^, for all t > tu (and t < —tn 
in the big crunch branch), where tn is the Hagedorn time. 

The case 6^ > 1 

We cannot use any more a function of the form (4.7) to map the Friedman equation 
of the present case to the pure de Sitter one. We thus actually solve eq. (2.5) in a more 
straightforward way. In the gauge = 1, one can immediately express a monotically 
increasing time t as a function of a : 

, , r udu 1 / 2Aa2-l \ 

tia) = / — - + ti = — z= arctanh — — — : , a > , 

Jo VAn4 _ ^,2 + ^2/4^ \2^/Xy/Xa*-a^ + 5^/4:\ J 

(4.29) 

where the argument of the square root is never vanishing and we have chosen the integration 
constant ti = — (4A)~^/^arctanh5^^. Inverting the function t(a), one finds 



a{t) = ^= \ 1 + J 51 - 1 sinh(2v^t) , t > U = y= arcsinh ^ . (4.30) 

v2A V V 2vA y5| 



The cosmological evolution described by this solution starts with a big bang at t = ti, while 
for large positive time, it is inflationary. It thus looks similar to the case 6^ < 1, both at the 
beginning of time and for large t, when the first order transition has occurred. However, in 
the present case, the evolution from the small to very large scale factors is smooth, (see fig. 
6). Actually, the solution (4.30) has an inflection point arising at t = ti^t, where a{tint) = ctinf, 



t... = ^arcsinhy|^, «... = ^ | ' (4-31) 

associated to a second order phase transition. As in the case 6^ < 1, the solution should be 
trusted as soon as the temperature is below Th- 

For completeness, we signal that another solution to eq. (2.5) is found under the time 
reversal t —t and is thus monotically decreasing. A de Sitter like universe is contracting 
since an infinitely past time and evolves in a smooth way around t = ti^f toward a branch 
that ends with a big crunch. 

The case 6^ = 1 
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Figure 6; Cosmological scenario of an expanding RW universe in presence of thermal and non trivial time- 
dependent moduli effects above some critical value. A phase of expansion starts from a big bang and evolves 
smoothly toward an inflationary one. This describes a second order transition. 



Taking the hmit — >^ 1 in the solutions found in the cases 5^ < 1 and 5'^ > 1 reaches the 
same static solution, 

a{t) = ao where ao = . (4.32) 

V 2A 

It corresponds to a universe with an S'^-space of constant radius. Actually, the 6^ — * 1_ limit 
amounts to send the big bang to an infinitely past time (see eqs. (4.14), (4.17)), while keeping 
the instant the transition occurs at t = 0. Similarly, the 6^ 1+ limit sends ti —oo (see 
eq. (4.30)), while keeping ti^f finite. In both cases, the solutions are thus "stretched" and 
converge toward a constant. In the two cases, it is however possible to redefine the big bang 
time at a conventional t = by a 5|.-dependent shift in the definition of time. Taking only 
then the limit 5^ — > 1 reaches solutions where the dynamical part of the solutions are not 
sent to infinity any more. 

To find these new solutions, one can write the Friedman equation (2.5) in the form 

{adf = A(a2 - aof , (4.33) 
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and derive two monotically increasing solutions for the time as a function of a, 

t{a) = -j= / , < a < ao , 4.34 

and 

1 r udu 



t{a) = — 2 ' « > oo . (4.35) 

Integrating explicitly these expressions and inverting the functions t{a) reaches then the two 
cosmological evolutions : 



a(t) =aoVl-e-2v^* , t>0, (4.36) 

and 

a{t) = aoVl + e2v^* , teR. (4.37) 

The first one starts with a big bang arising at t = 0, converges quickly toward the static 
solution, (and should be trusted as soon as T(t) < Th). The second is a deviation from the 
static solution toward the inflationary phase. Two other monotically decreasing solutions 
are found under the time reversal t — > — t in eqs. (4.36) and (4.37). The four solutions are 
drawn on fig. 7. 



5 De Sitter deformation by moduli fields 



In this section, we generalize the thermally deformed de Sitter solution by including the 
effect of the non trivial time dependence of some moduli fields. To be specific, our choice of 
matter action to be added to (2.3) is 



Ml , , 



(5.1) 



where the Xi's are the moduli fields with neither potential, nor interactions with the thermal 
system considered before. Assuming the RW metric (2.1) and the presence of the thermal 
corrections discussed before, the modified MSS action becomes 



iff 



dt Na-^ 



_1_ fd 

'm [a 



-i-A- 



4Aa4 6iV2 



a^d 



(5.2) 



where S'^ is a new positive constant. The energy- momentum tensor can be expressed in 
terms of 



1 



P 



36'^ 



Un^a^N^ V 4A 2^' 



(5.3) 
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Figure 7; Cosmological scenario of an expanding RW universe in presence of thermal and non trivial 
time- dependent moduli effects at a critical value. A phase of expansion starts from a big bang and converges 
quickly toward a static universe. After a very long time, the universe can either die in a big crunch or enter 
into an inflationary phase. 



^= hr-x + oxn • (5-4) 



127rV4iV2 \AX 2 
The equation of motion of Xi is trivially solved, 

X.^^^. (5.5) 

where 2^/26^. /3\ is an integration constant, whose form is chosen for later convenience. 
Defining 

Sl = j:^l.. (5.6) 

i 

the moduli contributions of p and P become 



p 1 



To solve the Friedman equation (2.5) in presence of the extra 1/a^-term induced by the 
motion of the Xi fields, we follow the method implemented in the previous sections. We look 
for a real constant n' in the definition of N{t), 

1 + A /a^[t) 
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so that the extra contribution is cancelled. This can be achieved by choosing n' such 
that 



V{n')^n'' + n- + ^-ln'-^^5l = Q. (5.9) 

This polynomial always has a real positive root (see the appendix). Using it, eq. (2.5) is 
taking the form 

'aV 1 + k' k'^ + k' + 5'^/A 



a) ■ <"°' 

The above equation has the same form with that of the 1/a^-thermally deformed de Sitter 
one discussed previously, up to a shift in both coefficients of the 1 / a^-curvature term and the 



1/a^-thermal term. Using a simple rescaling of the radius a (by a factor Vl + k'), one can 
immediately deduce the solutions from the previous section in terms of the time variable t. 
This implies in particular that the analogue of the order parameter 6^ of the pure thermal 
case is now the positive quantity 

{1 + K'y 27 «:'(! + K')' ' 

where these two definitions of A are equivalent as a consequence of eq. (5.9). In the appendix, 
we shown that the condition A < 1 is equivalent to the statement that the parameters 6!^ 
and 6^ satisfy the inequalities, 

6'S < 1 and 6l < h{6'S) , (5.12) 
where ^ 

This defines a domain in the (5^, 5^)-plane shown in fig. 8. At the origin of this domain, 
one has A = 0, while A = 1 on its boundary. Outside, this region, on has A > 1. 

The case A < 1 

Eq. (5.10) admits the cosmological evolution 



a{i) = M'\Je' + cosh^(v/At) , (5.14) 

where 
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Fi gure 8: Phase diagram in the (5!^ ,6'^) -plane. In the domain (5^ < h{6!^) defined in eq. (5.12), a 
first order phase transition between two cosmological branches can occur. Outside the domain, there is a 
solution describing a second order phase transition. On the critical curve (5^ — h{5'^), there are cosmological 
evolutions asymtotic to a static solution. Furthermore, one can introduce the temperatures < T- in the 
interior of the domain to replace the two parameters (5^ and (5y . On the boundary curve, these temperatures 
converge and their common value can be considered as a function Tct(5t) Tcxi^V) ■ 



To proceed, we switch to the more intuitive N = 1 gauge. This requests a redefinition of 
the time variable dictated by eq. (3.9) : 



t = I dv\l — — — V V where e" = , (5.16) 







An exphcit form of t as a function of t can actually be determined, however we find more 
convenient to work with the above integral representation. The scale factor a{t{t)) takes a 
simple form in terms of the function t{t) found by inverting the definition (5.16), 



a{i{t)) = £' + cosh^ {y^i{t)) . (5.17) 
Some comments are in order : 

• The slope dt/dt is positive everywhere, and for large t one has t ~ t. This means that the 
moduli deformation of the thermal de Sitter solution appears as a mild redefinition of time 
and normalization factor. 

• An Euclidean solution is defined via the analytic continuation, t = —if i.e. t = —ir, 



a^(f (r)) = Ar\/^M^^J^JS(7Xf(^ with r = - ^ dvj ^"^'N^^O + _ .^^^gx 
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The range of Euclidean time is 

-7r/2 < v^f < I.e. n < r < , (5.19) 

where Tj is given by 



dvJ . (5.20) 

--^ \ cos^{VXv)+e' + e" 

The boundary terms in eq. (3.3) vanish. 

• The Euchdean solution admits, (as in the pure thermal case), two analytic continuations, 
one for each boundary. The one associated to r^, ^/Xf = — 7r/2 + i^/Xi i.e. VXt = 
VXTf + i^/Xt, gives rise to a sinh-type cosmological solution. 



i(i{t)) = M'Je' - sinh2(VAt(t)) with t = - f dvJ ^' ^^'^^'(^^) . (5.21 
^ Ji \ e' + e" - sinh^(VAw) 



The range of time is 



arcsinhvV < t < «.e. = - / dv^l — ^ < t < 



.J\ y_ arcsinWZ V + — sinh^( J 

(5.22) 

The full scenario is thus similar to the one presented in the previous section. There 
is a big bang at t = tj and the universe expands till t = 0, before either contracting or 
experimenting a first order transition to an inflationary phase. At the transition, the scale 
factor is jumping. 



Thus, the present moduli deformation is generalizing the pure thermal one discussed previ- 
ously, (see fig. 3). 

The case A > 1 

Eq. (5.10) admits the t-dependent solution 



a{t) = V^^T^ a/i + sinh(2v^^ , t > U = ^ arcsinh^^= , (5.24) 

V 2A 2v A V A — 1 

we want to express in terms of a time variable associated to the natural = 1 gauge. 
Integrating eq. (3.9), one can define t as a function of i by 



t = I dv 



sinh(2yAt;) + 1 ^5) 



^ v/A^ sinh(x/2Ay) + 1 + 1^ 
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where 



t > U 



dv , 



y/K^ sinh(2v^i;) + 1 

2k' 



\^ v^A^ sinh(v^AtO + 1 



and rewrite 



a{i{t)) 



2A 



VA - 1 sinh(2yAt(t)) . 



(5.26) 



(5.27) 



As in the previous case, the slope dt/dt is always positive and t ~ t for large t. The solution 
is thus a mild deformation of the pure thermal case. It starts with a big bang aXt = ti and is 
monotically increasing. For large time, the behavior is inflationary. The evolution can again 
be interpreted as a second order phase transition since there is a unique inflection point, (see 
fig. 6). To show this, one can express dt/da as a function of a from the Friedman equation 
(in the = 1 gauge) and take a derivative. One finds that the inflection point arises at the 
radius ai^f, where 



Xa.^f > satisfies x —x — —6,, = . 



51 
27 ^ 



(5.28) 



Finally, a monotically decreasing cosmological solution is found under the time reversal 
t -t. 

The case A = 1 

The static solution obtained in the limit A ^ 1 in the previous cases is 



a{t) = Co where ao 



1 + k' 
2A 



(5.29) 



Beside this constant radius universe, we again have the two t-dependent evolutions 

a{i) = aoVlTe^^yx*. 
To consider them in the = 1 gauge, we integrate eq. (3.9) in each case and find 



(5.30) 



i{t{t)) = aoVl- e-2v^*W where t 



dv. 



V + ^ 



(5.31) 



for t > i.e. t > 0, and 



a(i{t)) = ao 



Vl + e2v^*(*) where t = / dv. 

Jo A 



o2^/Xv 



1 + + e 



(5.32) 
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for arbitrary i and t. As in the cases A 7^ 1, one has di/dt > and this redefinition of time 
does not change the quahtative behavior of the solutions, (see fig. 7). They are monotically 
increasing as in the pure thermal case, while two other solutions obtained under t ^ —t are 
decreasing. 

5.1 Other parameterizations of the solutions when A < 1 

The temperatures at the 1^^ order transition 

We would like to replace the parameters 5^ and 5'^ that appear in the scale factor solution 
by the temperatures T-t at t = O-t. From the definition (4.6) (with 5^ replacing 5^) and the 
relation (4.3), using the fact that aT = a+T!^ = a-T_, one has 



where the temperature has been defined in eq. (4.21). On the critical curve that is 
delimiting the domain (5.12), one has 6"^ = h{5'^). This implies A = 1 and thus T+ = T„. 
It is shown in the appendix that in this equals defined in eq. (A. 4). One 

can then deduce the critical temperature T^t on the boundary of the domain (5.12) by 



ploted on fig. 9. The temperature on the critical curve can also be considered as a function 




(5.33) 



T,,{5'^) = T^{6'lh{6'^)) (see fig. 8), 




(5.34) 



of (5^. This can be done by expressing 5'^ = h ^(5^), 




(5.35) 



where 




(5.36) 



Defining = T±{h-\6l),6l) (see fig. 8), one obtains 



l-{iA + BY/^ + iA-By/''-l) 
((A + 5)1/3 + (A -5)1/3 +1)2 



2\ 1/4 



(5.37) 
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Figure 9: Critical temperature as a function of 5'^ . It is defined on the curve S?. ~ h{5'^). 



Since T+a+ = T^a-, one has 



\2 ' 



(5.38) 



{T+/T_+T_/T^ 

a relation that generahzes eq. (4.21) and the duahty ^ . Actually, for eqs. (5.38) and 
(4.21) to have precisely the same form, one can introduce a mean temperature = T+T_, 



T' =T 



5'S 



4/t'2 + 4k' + 6'^ 



1/4 



(5.39) 



where has been defined in eq. (4.21). However, this definition is less useful than its 
counterpart in the pure thermal case. This is due to the fact that depends on the point 
(5^, 6"^) of the parameter space, while in the pure thermal case, 5^ = = so that the mean 
temperature Tm is independent of the remaining parameter 6!^. 

We may express the evolution of the scale factor a{t{t)) in each branch with T±, 



and 



a{i{t)) 



a{i{t)) 



V2A(1-^) 



1 + t:^ ^ cosh(2v^t») , t > 



1 + -A — cosh{2Vxi{t)) , t < , 
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(5.40) 



(5.41) 



where k' has been expressed in terms of A that involves T± and T„ 



1 1-T^T?/T^ 

where A = 

1-A {T+/T^+T_/T^ 



l + n' = where A = ^. ~ ,Z 2 ■ (5-42) 



To obtain the above expression for A, we have used eq. (5.11) to express 5^ in terms of A 
and k'. Then we combine the resuh with eqs. (5.33) and (5.38). 

Under the duahty T+ T_, the two cosmological solutions a{t{t)), for t > and t < 
are interchanged : 

( T+ ^ T_) [a{i{t)) for t > ^ a{i{t)) for t < o) . (5.43) 

For completeness we display the constants A/"', e' and e" appearing in eqs. (5.17), (5.16) 
and (5.21), in terms of T±, 



^K\-A) V 1 + {T^IT-f ' (T_/T+)2 - 1 ' (r_/T+)2 - (T+/r_)2 • 

(5.44) 

Parameterization with 5'^ and 5'^ 

Let us determine explicit expressions of the scale factor in terms of 5'^ and 5^, when A < 1. 
This will be useful in the next section, when we study the probability transition. Since the 
parameters 5^ and 5^ are chosen in the domain (5.12), the polynomial V in eq. (5.45) has 
three real roots, (see appendix). Noting that 'P{—\a±) = 0, as can be checked with eq. 
(5.11), the polynomial V is actually, 

V(x) = {x- k'){x + Xa^_){x + Xal) . (5.45) 

Using Cardan's formulas, one finds 

k' = 1 (j4-36'Scos (l) - 1 ) , (5.46) 



3 V V ^ V 3 



e±27i 



3^.1-V4-35tCos(^— ) ) , (5.47) 



where 

1661 + 96'^ -8 



e = arccos I ^^^7,^^ 1 . (5.48) 

' (4-35^2)3/2 / 
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We note that all values of 6 in [0,7r] can be reached. At the origin (5^,5^) = (0,0) of the 
domain (5.12), one has 6 = n, while 6 vanishes on the critical curve 5^ = h{6'^). The constant 
A/"' of eqs. (5.17) and (5.21) can be determined by expressing y/1 — A from the ratio a^/ 



The relations = Af^e' and k' = XJ^f^e" give then 

_ 1 - cos (^) ^„ ^ V^^cos(f)-l ^ 



v^V^4^^Wsin (f) ' ysyi^^Wsind) 
For completeness, the temperatures at the transition are also found to be 



^5' 
2"r 



T± = T^.^= ' (5.51) 
- v/4^^cos(^) 



5.2 First order transition probability 

When A < 1, the transition probability at t = between the big bang branch and the de- 
formed de Sitter cosmology can be estimated by computing the Euclidean action. Changing 
of integration variable dr = daE/dE in the Euclidean analog of the action (5.2) gives, on 
shell. 



where V is given in eq. (5.45) and enters in the Euclidean equation of motion 

(5.53) 



^qe J a% 

we have used. Clearly, SEeff is negative. To explicitly evaluate it, we prefer considering its 
integral form on f and use eq. (2.5) in Euclidean time to be on shell. 



|-o ^ /I 5'"^ \ 



(5.54) 

\U,~E ^AU'e J 

2\/A 

In this expression, a£;(^) and NeIt) = dr/df are referring to eq. (5.18). The result can 
be expressed in terms of complete elliptic integrals of the first kind, K{k), and second kind. 
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E{k), 

^E.ff=-1^{'^-^ST) Wsm ^ \Eik) ^ r^^,_,__ ^ TTZ^K{k) 



(5.55) 



where 



, sin (I) 

To have a better intuition of the behavior of this action as a function of the parameters S!^ 

and (5^, we concentrate on the boundary of the domain (5.12) : 

• On the side 6"?. = corresponding to the pure thermal case, the action becomes 




Se., = ^ ' ""T^ ( E{k) -{l-\/l-S'^] K{k) ) , (5.57) 

where 

2(1 - 5'2)i/4 

k= ; ^' (5.5^ 



1 + 



To derive these expressions, one can use the fact that in eq. (5.46) vanishes and is thus 
giving rise to an identity satisfied by 6. Eq. (5.57) reproduces the resuh of [10]^. At the 
origin (5^,5^) = (0,0), one has k = 1 and £'(1) = 1, so that the resuh — 1/3A of the pure 
de Sitter case is recovered. At {6'S,Sl) = (1,0), k = and E{0) = K{0) = 7r/2, so that the 
action vanishes. This is consistent with eq. (5.52) since a+ = a_ and 5^ vanishes. 
• On the critical curve 5^ = h{5'^), the action is taking the form 



, ^1 + ^1 - f (^1 - |i? - 1/2^ 

If it is explicitly negative for 5^ < 1, it vanishes for 6'^ = 1 only. This shows that even if 
a+ = a_ on the critical curve, the term proportional to 5^ in the integrand of eq. (5.52) 
contributes, due to the fact that the denominator vanishes. The probability p refers in 
this case to the transition (after an infinite time) from the big bang branch (5.31) to the 
inflationary one (5.32). The quantity (1 — p) refers to the transition (after an infinite time) 
form the big bang branch (5.31) to its time reversal i.e. describing a big crunch. This is 
the case since considering the solutions (5.17) and (5.21) of the case A < 1, and performing 



''Up to minor misprint errors in [10]. 
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shifts on the origins of times, give rise to (5.31) and (5.32) in the hmit A 1_. 
• The side 5^ = is corresponding to the pure moduh deformation of the de Sitter solution. 
From eq. (5.59), one finds that at the corner (5^,5^) = (0,1) of the domain (5.12), the 
action is taking the value — 7r/9A. Since this result is close to the — 1/3A at the origin of the 
(5^, (5^)-plane, the dependence of the action on 6^ seems to be mild along this axis. 

This remark happens to remain true for arbitrary fixed 6'^. To understand this, one can 
note in eq. (5.54) that the S'^/a^ contribution of Xi in the action SEeg has disappeared on 
shell. Thus, the only dependence on 6"^ occurs through the mild deformation e" appearing 
in Ne- This is confirmed on the 3-dimensional plot of fig. 10, where SEeff- appears almost 
constant when 5^ varies from to h{6'^), at fixed 5^. 




Figure 10: Euclidean action SecB as a function of the radiation parameters'^ and the moduli deformation 
S^. The variable ((5^,(5^) spans the domain (5.12). The dependence on is weak, while the absolute value 
of SecB decreases when 8'^ is switched on. (We have chosen A = 1 on this plot.) 

From fig. 10, one can also see that the transition amplitude p oc e"^'^-^"-^ is larger and 
larger when 5^ decreases from /i~^((5^) to 0, at fixed 5^. However, the magnitude of p highly 
depends on A. For instance, p may not be large even for small 6!^, due to the 1 / A-dependance 
of Secs, for large A. 
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At fixed A and 5^, switching on 5^ increases e' (see eqs. (5.50) and (5.48)), and thus 



2\ti\ (see eq. (5.22)), the range of time along the big bang / big crunch evolution. From 



increasing 5^ makes this cosmological branch more and more viable. 

Similarly, at fixed A and 5^ 7^ 0, switching on 5'^ has the effect to increase 2|tj|. However, 
T_ varies from to T^y. , with the risk for the universe to be too hot. The case = has a 
different behavior. When 5'^ grows up from to 1, 2\ti\ increases as well but the temperature 
T_ varies from +00 to Tm-^ Along this axis, 5^ is thus making the big bang / big crunch 
branch more viable. 

Finally, when we reach the limit case corresponding to (5^,5^) sitting on the critical 
curve, after the universe is born with a big bang, it converges quickly toward a static state 
for a very long time, i.e. an 5''^ of constant radius 



where oq has already been defined in eq. (5.29). However, it finally enters in an inflationary 
phase or dies in a big crunch. Also, one can note that in the case A > 1, the big bang 
evolution is longer and longer when A 1+, since ti —00, (see eqs. (5.24) and (5.26)), 
but the final state on the universe is always inflationary. 

6 Conclusions 

In this work, we analyze some aspects concerning inflationary solutions and possible tran- 
sitions between different cosmological behaviors. Our starting point is a bare de Sitter 
background deformed by the presence of a thermal bath and the motion of moduli fields. 
The thermal bath summarizes the effects of light degrees of freedom in a consistent way and 
shows up as a Cnjo?' correction in the MSS-action. The moduli contribution gives a Cm/cl^ 
correction to the energy density. 

A systematic analysis of the system involves an almost triangular domain in the defor- 

*Note that r_ +00 when (5^ = and (5^ ^ 0, while r_ = when (5^ = and 7^ 0. This ambiguity 
in the definition of in the pure de Sitter universe is due to the fact that taking (5^ = (i.e. k' = 0) or 
(5^ = are operations that do not commute (see eq. (5.39)). 



eq. (5.51), larger and larger 5^ 



imply also lower and lower temperature T_. As a result. 




(5.60) 
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mation parameter space (5^,5^). We solve the gravitational equations in all cases and find 
that a first order phase transition can occur inside this domain, while a second order one 
arises outside. 

In the second case, a smooth transition always occurs between the big bang cosmological 
evolution (as soon as T(t) < Th) and the inflationary behavior. In the first case, we calculate 
the transition probability between the big bang cosmological evolution (where T_ < T{t) < 
Th) and the inflationary branch (where T{t) < r+). The two solutions are connected into 
each other via a gravitational instanton allowing a double analytic continuation. We find, 
the existence of a temperature duality T+ ^ T_ that interchanges the two solutions. Inside 
the domain of the parameter space, the origin (5^,5^) ^ (0,0) corresponds to the pure de 
Sitter case where the big-bang branch disappears and the instanton connects the inflationary 
universe to "nothing". In the two other extreme cases, i) the pure thermal deformation 
(5^, 5^) ~ (1, 0) and ii) the pure moduli deformation (5^, 6^) ~ (0, 1), the big bang branch 
remains almost static for a very long period of time. However, the probability transition is 
minimal in case i) and maximal in case ii). Thus, the late future of these evolutions is most 
probably contracting (till a big crunch occurs) in case i) , and inflationary in case ii) . 

The above analysis and results can be applied to more complex systems where the radia- 
tive and temperature corrections are effectively taking the form of cosmological, curvature 
and radiation terms during the time evolution [12]. These richer models occur in string effec- 
tive no-scale supergravities i.e. associated to = 1 string compactifications where = 1 
supersymmetry is spontaneously broken. First and second order phase transitions are again 
involved [12]. The stringy origin of these models is giving us the possibility to go beyond the 
field theory approximation. In particular, it is possible to study the big bang cosmological 
evolution above the Hagedorn temperature [18] as well as "the stringy analog of the wave 
function of the universe" [19]. 
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Appendix 

We would like to find in this appendix the phase diagram in the ((5^, 5^) -plane that is 
delimiting when the solutions (5.17) and (5.21) for the scale factor are valid, and when it is 
instead the evolution (5.24) that is relevant. The two first solutions exist when A < 1, while 
the last one occurs when A > 1. 

First of all, we note that the polynomial V defined in eq. (5.9) has three roots, whose 
product is 45^/27. If two of them are complex conjugate, the third one is thus real positive. 
If on the contrary they are all real, one or three of them must be positive. However, since 
the sum of these real roots is —1, only one is positive. In any case, we always have a single 
real positive root, k', and eventually two other real negative ones. 

For 5'^ > 4/3, the derivative V of the polynomial defined in eq. (5.9) is everywhere 
positive, so that V has a single root. When 

S?<1, (A.l) 

the derivative V admits two roots, 

= (^1 T /l^j ■ (A.2) 

One always has V{k'^) < 0, while V{k,'_) > (so that V has 3 real roots) if and only 
if 6^ < /i(5y), where h{6''^) is defined in eq. (5.13). Actually, this condition implies in 
particular that h{5'^) > 0, which is equivalent to S'^ < 1. As a conclusion, the polynomial V 
has three real roots, (one positive and two negative), if the inequalities (5.12) are satisfied, 
while it has only one real root in all other cases. 

Let us determine when A < 1. This condition can be translated into a degree two 
inequality for n' that admits solutions if the condition (A.l) is satisfied. In that case, A < 1 
is equivalent to having 

^/\— — ^ — ' (A*3) 
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where 

= ^ (±y/4-35^2 - . (A.4) 

Since having < is equivalent to satisfying 5^ > 1, the single real root k' that V has 
in this case is outside the range /^a+I' since it is positive. Thus, 5'^ > 1 implies A > 1. 

Let us then concentrate on the case 5'^ < 1. We note that when the parameters 5^ and 
(5^ are chosen on the critical curve defined by 5^ = one has 'P(/tA+) = 0. This means 

that the single real positive root n' of V is then precisely i^aX case, we have A = 1. 

If instead we take 5^ > h{5'^), the only real root of V is positive and satisfies hi' > 
due to the fact that V'{x) > for any x > 0. We thus have A > 1. On the contrary, if 
5^ < h{6'^), the real root of V which is positive satisfies k' < k,'^_^ for the same reason, so 
that we have A < 1. 

As a conclusion, the solutions (5.17) and (5.21) corresponding to the case A < 1 arise 
when the point (5^,5^) belongs to the interior of the domain (5.12) of fig. 8. Outside this 
domain, the evolution (5.24) associated to A > 1 is the correct one. On the critical curve 
(52 = h{6'S), one has A = 1. 
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